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PREFACE 


The purpose of this research was to examine the instantaneous 
divergence of the rradient wind at individual points in circular and 
sinusoidal isobaric patternse 

This work was conducted at the U. S. Naval Postgraduate School, 
Wonterey, California, during the period December 1949 to May 1950. 

I wish to express my appreciation to Professor William De Duthie 
and to Associate Professor A. Be Mewoorn for their advice and ruidance 
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DRODUCT ION 


The variation of pressure with time in any meteorological situation 
has long been recognized as one of the basic foundations upon which pro= 
gnosis and forecast must rely. As a result, a considerable amount of 
meteorolozical literature has been devoted to the tendency equation in 
its various formse Both the advection term and the diverrence term have 
been extensively investigated from the theoretical standpoint. Potterssen 
[9]in a recent paper devised a simple method for quantitative ovaluaticn 
of the advection term. As yet no such method has been devised for the 
diversence term. Bjerkmes and Holmooe (4] considered mass diverrence 
instead of investiratins the two terms separately (as is dono here). 

Their approach was dynamic in nature, beins associated with the equations 
of motion and isobars or isohypses on unper-level maps. 

‘ost investigations of this subject have considered the averase dis~ 
tribution of divergence over an area. A study of the instantaneous velccity 
diversence in circular and sinusoidal patterns is considered hore, and points 
as well as areas of convergence and diversence comparede 

The first chapter deals with a derivation and a brief qualitative 
discussion of mass divergence. The mmericeal example, however, considers 
density a constant, thereby siving tho nroduct of density times the diver- 
gence of the gradient wind. Throughout the followins, unless mass diverrence 
is specifically specified, velocity diversence is intended. Likewise, only 


horizontal divergence is being considered in this study. 
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Worizontal divergence of any vector-field is a scalar function of 
position and time. At a fixed time therofore, horizontal divergence is 
a function of position alone. Thus, if the vector=field Yy is known at 
time (t) , say at tho time of the woather map, tho divergence may be deter-= 
mined from the vector=field at that timee The horizontal divergence in 
this investigation is broken dovm into two parts: (1) the wind-~sneed 
times the diverrence of the unitevector tanrent to the horizontal strean- 
lines; and (2) the scalar product of the wnitetangent vector and the 
gradient of the horizontal wind-speed. Five different isobar patterns 
are investigated: (1) stationary circular low with constant pressure 
eradient; (2) circular low, with constant pressure sradient, moving east- 
ward 10 lmots; (3) circular low, with constant pressure gradient, moving 
eastward 30 kmots; (4) stationary circular low with pressure cradient a 
function of radius of curvature of isobarss; and £52 Sinusoidal wavese 
Horizontal diverrence is computed at various points in each of the patterns 
and the relative magnitudes notede In all the circular cases investigated, 
there is convergence to the east of the meridian running throurh the center 
and divergence to the west. The areas of maximum converrence and diverrence 
are on the innermost isobars along the latitude circles throurh the centere 
The areas of instantaneous convergence and divergence are qualitatively the 
same as the areas of average convergence and diverrgonce as stated by Bierkmes 
and Holmboe C4] . The general areas of maximum and minimum convergence end 
divergence in circular lows are in qualitative agreement with the findines 


of Petterssen [8 , ppe 236-237]. 
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. theorotical example is computed using a circular low at tho 
surface with various average sinusoidal waves alofte The changes in the 
contributions to the surface pressure tendoncy with che differcnt wave 
pattorns aloft are recordod in Table 5. 

An invostication of the circular systems in this study shows that 
the isovels bear a close relationship to the critical eccontricity statod 
by Bjerlmes and Holmboe [4] for circular systems. This would be expected 
however, sinco the isovels are lines of constant ( aA +VK) wie is 
vroporvional to tho transnort capacity of isobaric channels as stated by 
Bjerlmes and Holmboe [4] .« 

A computation from the 6850, 700, 500 and 300 millibar charts of 15002, 
19 October, 1949 is made and tho resulting contributions to the surface 
pressure tendency are compared with the observed tendencies. The order of 
magnitude of the contribution of the divergence term from the layers ine 
vestigated is in agreement with the results of Baum [2] . 

Gradient flow is assumed throushout. The wind of course is not always 
cradiont in nature even in the free atmosphere, although the latest investi- 
gations have shown that the difference vetween the sradient vrind and the 
observed winds is usually very small at about 1000 meters above the srounde 
yinen the pressure gradient is changing rapidly, however, the relationship 
between the observed wind ond the pressure gradient is considerably modified. 
This is due to the fact that the motion is then not under balanced forces. 
Under tnese conditions there is a velocity component along the isallobaric 
sradicnte A computation of divergence assuming gradient flow in this case 


would generally be in error. 
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The synoptic evidences of convergence and divergence are discussed. 
In general it can be assumed that widespread cloudiness implies widespread 


converrence and clear skies imply divergenco over the area. 


(4) 
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specific momentum in gsradient flow as the sealar product of a vector parallel 


Ie THE MASS DIVERGCUCE EQUATION 


Theoretical derivatione 


Let the winc be expressed in vector notation as 
Y= V5 — (us) v (1.1) 
Multiplying equation (1.1) by CA and collecting terms fives 
(A+vk)ey = Aeu,. = -vPxk (aa 
Solving equation (1.2) for CY we cet 


—/ 
C@v= ~ VPXK (A + vk) (153) 
Talking the horizontal divergence of both sides gives 


kl = 
VeCWe =VPKik* V(t Vm) + (ave) V-VP RK 


which reduces to 
-| 
Vev = -VPXK + Viaeve) (164) 


since V-VPXK = Oe 


Equation (1e4) therefore gives the mass divergence or divergence of 


to the isobars and the gradient of the quantity ( At+VK), 


(5) 
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Expansion of the right hand side of equation (1.4) fives 


=e, 
VCs UAtvel VPXKK o> 9 (Cheval) (145) 
i -2 
=(A+ve)TPXK*VA + (Atv) VPXK: Vive) (146) 


oP huts - Valen, VPKE ‘AVY + foe. Ew, 


(aA +Vvk)* (A + VK 


iA 


Since VAis along a meridian, it may be written as VA= oy¥Z . 
Therefore, the first term on the right side of equation (le7) can be 


further simplified. Making use of the fact that K= % © and 
dd + 
dy = fy (kw Sin Se} Pa 
2 Je 
VA reduces to 


VA = Act od 


Rouation (le7) may therefore be written as 


=a —2 . ~2 
Vs CY = (4 thew OM AtVK) AVPKK + Z + (AtVK)VPXK KV 


~2 
Oh PV YP ie. Vin” © (1.8) 
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2e Qualitative Interpretation, 

A qualitative oxamination of equation (1.6) shows that for straight 
parallel isobars running east west, the total mass divergence is exactly 
zeroe Tho first term on the right hand side is zero because the two vectors 
are perpendicular; the second, because kK is zero for straight isobars, If 
the isobars have any northward component, there is convergence, and any 
southward component, there is divergence. This convergence or divergence 
reaches a maximum when the isobars run duo north or south, other things 
being equale 

Consider next the "Cartwheel", with vk caual to a constant. From 
equation (1.6), it can be seen that the complete divergence is given by 
the tern (AsVK) UPXK -VA . Sinee YA= Weai Ce Od , the varietion 
in divergence in middle latitudes for this case is largely determined by the 
direction and marcnitude of the vector VPA K, ~or the variations in A, @ and 
Cot @ are very small, ‘The area of convergence is to the east of the merid- 
ian throush the center, and the area of divergence to the weste The magnitudes 
at points equatorward of the center are slightly larcer than the values for 
corresponding points on the poleward side of the latitude circle through the 
center 

lor a stationary circular low with constant pressure gradient, the 
numerator of the first term on the right side of equation (1.7) is zero on 
the meridian running throuzh the center of the low, and reaches a maximum 
along the latitude circle through the center of the low if we neglect the 
very small variation of VA in middle latitudes. It is negative to the 


eastward of the meridian throuzh the center and positive to the westward 


(7) 
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of the meridian. With density also a constant, the second term of (le7), 
as will be seen from equation (1.10), is zero on tho moridian through the 
center, reaches a maximum on tho innermost isobar closo to the latitude 
circle through the center, is nesative to the eastward of the meridian 
throush the center and positive to the westward of the meridian. 
3e $A Vumerical Example. 

For a numerical investigation of the stationary circular low with 
constant pressure gradient, isolines of the quantity (aA +VK) rust be 


dravm. This is done in the following manner. Set 


=| 
(A+tvK) = CONSTRAT = W 
oe 5 — aT (led 


Substituting this value of vk in the gradient wind equation gives 
oP 
AV +V(R-A) =¢ am 
Pome (28 (1-10) 


€ IN ° 


With constant pressure gradient and the variation of density in the hori- 
zontal neglected, the isolines of (A a are isovels, and are labeled 
with the appropriate value of the wind speed in Figure 1. Choosing an 
isovel to be plotted (10 meters per second for example), the corresponding 
value of VW can be computed from equation (1.10). With v end Vv lmom, 
assigning the values of A that occur in the low will determine the correct 
values of Ki « These A and K: values can then be plotted on the system 


using latitude circles and radii of curvature of isobars as coordinates. 


Connecting this series of points will give the isovel desirede The isovels 


(8) 
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constructed for the stationary circular case are shown in Fisure le from 
Ficure 1 it is seen that there is convergence to the eastward and diverzcence 
to the westward of the meridian through the center, end the maximum areas of 
convergence and divergence are on the innermost isobars alongs the latitude 


circles through the center. Table 1 sives the relative values of divergence 


computed at the positions indicated in Figure l. 


Value of 


Position Divergence x 10? 
A -e dl 
B ~e22 
C ~e52 
D ~o14 
E 015 
F 022 
G 039 
H 008 
TABLE 1, 


Since density is considered constant in this computation, equation (144) 
no lonser fives mass divergence but the product of density times the diver- 
gence of the cradient wind. 

Figure 1 shows an interesting comparison with the theory of Bjerkmes 
and Holmboe [4] on the critical eccentricity necessary for a scvationary 
low. Their results showed that a stationary low with exactly the critical 
eccentricity for all isobaric channels would have almost concentric isobars 
near the center and increasingly eccentric isobars toward the peripherye A 
glance at Figure 1 shows that the isovels have almost exactly this distri- 


butione We would expect them to have this distribution, however, since the 


(9) 
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Figure le 
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ail, 


isovels are lines of constant (A+VK), which, according to Bijerimes 
and Holmboe ([(4),are proportional to the transport capacity in gradient 
flow. Consequently, the isovels in Figures 1, 5 and 6 sive the occentri- 


city of the isolines of equal transport capacity. 
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Il. THS VELOCITY DIVEPGSHGl MQUATION 


1. Derivation of Formula. 
A glance at eavation (1.10) shows that unless the pressure sradient 
a! 
is constant, isolines of the quantit:: (A+V«K) are not isovels. In order 
to determine the diverrence from the isovels, the followinr notations were 


adopted. Set the wind equal to scalar y times vector % where t is 


defined as a unit vector varallel to the wind. 


Y=aVe (2.1) 


\ 


t =- VE XK (a2 


the horizontal divergence of the wind may be written, oy a simple 


vector=identity, as 
We Mice N/a to (23) 


Equation (203) breaks the horizontal divergence dowm into two parts: 
(1) the wind speed times the diverrence of the unitevector tanzent to the 
streamlines; and (2) the scalar nroduct of the unit-vector tanrent to the 
streamlines and the sradient of the horizontal wind speed. 

Investigating the first term on the richt side of equation (2.3), we 


find by substitutins the value of t given in equation (22) 


Va Mae? a Se Mel g Fo Kk 


= NOV PX K+ (-V PAK Vive ) 
(204) 


V4 


-V VPKK + Yj 


since V- VPX K=0, 
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Equation (264) may be further simplified by expanding the right 
hand terme 


VV-me-WVPik- Vie = VVPX - am 
lye 


wy VEX. 
va | vP/ 


=-WV- Vive (205) 
| VP/ 


2e Qualitative Interpretation. 

Equation (2.5) shows that with a constant pressure gradient, the 
term VV-t=oe A closer examination shows that for parallel isobars or 
for concentric circular isobars (not necessarily equally spaced) the term 
VV-& is still zero, since the gradient of {VP} is normal to the wind 
vector. For these three interesting cases then, the complete divergence 
of the wind is given by the term ti VV « 

VV-& essentially measures the spreading or closing of the isobars 
as can be seen from equation (2.5). ‘ith isobars spreading dowmrind, the 
term VV-% is positive, giving a negative contribution to the pressure 
tendency, This is in agreement with Scherhag's divergence theoren as 
stated by Baum [ 2] , that a horizontal flow divergence aloft is mere 
frequently associated with pressure falls than with pressure rises belowWe 

Examination of the second term of equation (23) shows that it is also 
related to the spreading or closing of the isobars. In general it is of 
onnosite sisn from the termVV-+% . Therefore, the two terms tend to cancel 


each other, and the total contribution to the pressure tendency is the residue. 


(13) 
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Neither term by itself then gives a completo picture of the divergence 
unless the other term is zeroe VY¥-% is zero for Gp constant, parallel 
isobars or concentric circular isobarse t%.VV is zero when the gradient 
of wind speed is normal to the wind direction. Im these rostricted cases, 
the complete divergence is givon by the appropriate term in equation (2.3). 
5e <A Numerical Example. 

A svationary circular low with the pressure gradient a function of radius 
of curvature of the isobars is constructed. 

The gradient wind is a function of four variables: A, @, Ri, and VP. 
Density has been considered constant in the horizontal throuchout this inves- 


tigation. Therefore, the wind can be written as 


V = F(A,R:, 28, ) (2.6) 


Since the low under consideration here is constructed with the pressure 
Gradient a function of radius of curvature of the isobars, equation (2.6) 


reduces to 
V=aF (A R:) (207) 


For plotting isovels therefore, choosing the value of vy for the desired 
isoline and assiming values to R:, determines the corresponding values 
for A. With values of A and R: determined, the isovel is dravm by 
connectinz the points by a smooth curvee Fisure 2 shows the results of 
the computations. Figure 3 is the pressure profile of the low under con~ 


Sideratione 
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4e Conclusions from the Variable Prossure Gradient Case, 

Firure 2 shows that the areas of converrence and diverrence are the 
seme as in the previous case, with convergence to the east and divergence 
to the wost. The relative mapnitudes changed, however. The computed 
values of the convergence to the east of the system were larrer than those 
computed in the previous case. Table 2 fives the relative values of diver- 


gence compuved at the positions indicated in Figure 2. 


Value of 
Position Divergence x 10? 


= 408 
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IIIe NUMERICAL INVESTIGATIONS 


le Movins Systems. 

Isovels aro arain plotted to detormine quantitative values of diverrence 
ina lowe A circular low with constant pressure fradient is constructedse 
Fastward movements of 10 and 30 Imots are considerede 

Since the K appearing in equation (le4) is the curvature of the tra- 
jectory of the air particle, it is necessary to relate it to the curvature 
of the isobars for plotting purposes. If the X-axis is chosen along the 
path of the center, and ansle y is the angle between the K-axis and the 


wind direction, then from Petterssen [8] the desired relationship is 


K= ke (1-€ Cs J (301) 





Firure 4e¢ 


(18 ) 





To simplify plotting, the following mothod was derived. Take an 
arbitrary point P.on an arbitrary isobar (Figure 4.). The line ccnnecting 
P and the center of the isobar is Ri. Drawa line from the center of the 
isobar in the direction of movement of the system. Construct a tanrent to 
tho isobar at point P and extend the tangent line until it intersects tho 
line in the direction of movement of the system. The angle between theso 
two lines is then, by definition, angle yp. Let af , measured in meters, 
represent the distance that point P is north or south of the center of the 


isobars. Then, by trigonometry, 


AY = Ri Sin OP = “Ree Sin Gee 


i 


I 


oy Gin (wv -fe) 


_ 
Eee 


{! 


ES Se — Ks AY (302) 


Substitutins the value of CoS Y above in equation (361) and solving gives 


Ria VS VV ea (363) 
2VK 
In order to plot the isovels, a value is assumed (10 MPS for example) 
which determines a value of NW in equation (1.10). Then choosing the values 
of A that appear in the system, the corresponding values of K could be de- 
termined. The correct values of A and Ri could then be computed from 


equation (303). 


(19) 
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2e Conclusions from Moving Systoms. 

Fisures 5 and 6 represent oastward moving systems with speeds of 10 
and 30 knots respectively. Increasing the eastward movement of the system 
increases the sradient of VY in the southern portion of the systom and de-= 
creases the gradient of wv in the northern portion. Considering the changes 
of the vector field YW, this is what we would expect. Since the isobar 
patterns were identical, the directions in the vector field were not chanzede 
The magnitudes, however, did undergo a chansee In the southern portion of 
the systom, the movement of the system added to the magnitude of the wind 
vector since it was in the same direction. In the northern portion, on the 
other hand, the movement of the system was opposite to the wind vector so 
deducted from ite The area of convergence is still to the east of the 
meridian through the center of the low and divergence to the west of the 
meridians 

The alsebraic sim of equation (le4) or (203) can be determined from a 
slance at Figures 1, 5 and 6.4 There the isovels cross isobars toward hisher 
pressure, there is convergence and where they so toward lower pressure, there 
is divergencee Table 4 gives the relative values of divergence computed at 
the positions indicated in Figure 5- 

Value of 

Position Divergence x 10? 

=28 
=e10 
=097 
~eO] 
e2l 
026 


029 
007 


qQ WMO Qt Se 


TABLE 3. 


(20) 
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de Investigation of Sinusoidal Waves. 

Four sine waves are dravm displaced 180 nautical miles at the trough 
and wedso lines (Plate Ie). Their wave lonzth is 1800 nautical miles and 
their amplitude 660 nautical milos. Let the%-axis point northward and 
the X -axis point eastward. Then the equation of the sine weve may bo 


civen as 


Va = Bi S/neex 





(205) 
The curvature at any point on the sine wave is given by 
a 
K = 4% 
Zin * 
a a 3 
ie ak 
a e 
= fa B Sin BX 
= ae 
[1 + (£8 cos 8x) | ‘ | (3.6) 


the value of K is computed at 10° intervals along the east-west axise 
Measuring the pressure pfradient, and computine the value of A at each of 
these points, the value of the rradient wind can be obtained fron equation 


(307) for cyclonic curvature and equation (3.8) for anticyclonic curvature. 
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VezaAR(I-Vi- #& ) (348) 


AER 
Using the values of tho wind computed at the various points, isovels are 


dravm as shovm in Plate I, 
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4.  offparigon of the Tersist -vv ani V¥-t . 
tho elretraic gicr of tho term +4 -Vv in weaeeés sal troushe with 


sinusoidal isobars is eiven ty Tirure 7. 


Fieure Je 


} « 


Lealine at the sien of VV-t in Ticvre u, it is evide>t &hat ‘the to 


isn 


verms tend to cancel each other since t\ryv are of opposite siens and are 


@ 
a 


both related to the pressure fradien%. 





The values of the two terms aro computed at the points A, 3, 7%, D 


and 5 on Plate Ie VV-Gis computed by the identity 


=A 
\/ te Niles) (lod ek D (209) 


as sucgestod by Beers [3], where L:7-% is siven by 


~2 
LiVet = 81 (RF) Sin (27%) Cos (27%) 


= (3% 10) 
- : 
[1 ba 70 (24 ) Ces* (2m K) 1% 


As can be seen from Table 4, t Vv is the dominant term in this parti- 


cular case, but both terms are of the same order of marnitudese 


Point t .VV X 10° VV-b xt 
A 1034 -1.07 
B 042 - 16 
C -126 094 
D 1633 = .% 
L = oD 61 
TABLE 4.6 
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IVe S¥VOPTIC ASPTCTS 


1. Avnlication to the Tendency Saquatione 
Petterssen [9, ppe 2-3], using the hydrostatic equation and the 


equation of continuity, derived the followine for the tendency equation: 


Ce - OF e PRRY wed Cy, (4.1) 


which for sea level roducos to 


“sh eV +P VY. (4¥2) 


A quantitative evaluation of the first term on the risht side of equation 
(4.1) was shovm by Petterssen [9] to be fairly simple. Vector YW is the 
isovaric=mean wind. It is obtained by observing the winds at equidistant 
pressure intervals (each isobaric layer containing the same amount of mass), 
surming, and dividing by the number of pressure intervals. Tho scalar pro- 
duct of this mean wind and the surface pressure gradient is then computed, 
giving the contribution of the transport term to the surface pressure ten- 
dency. The second, or divergence term, may be evaluated in a similar 


mannere Since 


pv-Y =RL V Ye (403) 


ate 
kt 


evaluation of the rirht side of equation (4.3) gives the contribution of 


the divergence term to the surface pressure tendencye 


(26) 
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2° Theoretical Investiration of a Surfaco Low with a Sinusoidal 
Favvern Above es 


Consider a surface low, with an averare cross=<isobar flow of TS? 
averaco wind speed (to the ton of the friction-layer) of 30 Imots, radius 
of curvature of isobar 180 nautical miles, and a friction=-laycr 100 mb thick 
Thon the contribution of the frictional inflow to the surface pressuro ten- 
doncy is approximately 12 millitars per 5 hourse 

Next consider respectively the points A, B, C, D and E of Plate I as 
being the averase diverrence nattern for a thickmess of 500 millibars above 
& point on the low being considered. Then Table 5 fives the total contri-g 


bution of the two layers to the surface pressure tendencye 


Millibars per 3 Hours 


Fr@tion 
Point Layer L-VV + vy¥-t% Total 
A 12 o1365 ~14.5 
B al H1420 ~2 0 
C 12 LT ec 29 2 
D 12 =-21,0 920 
B 2 13520 2540 


TABLE 5e 


(27) 





36 Computation of &-VV and VV-°t from Actua) Chartse 

The surface, 850, 700, 500 and 300 millibar charts of 15002, 
Octobor 19, 1949 are used. Isovels are dravm on all tho upper air charts 
as accurately as possible considering the scarcity of reports. The mean 
value of &-VV for the layers is then computed. Likewise, the mean value 
of VV°& is obtained. The contribution of the two terms to the surface 
pressure tendency is computed at six stations. Only stations with steady 


rising or falling tendencies are selected. The results are tabulated in 


Table 6. 
mb per 3 Hours 
Station t-VV VVet Total pt 
476 27 =45 ~18 ~19 
469 47 =44 3 1e4 
465 28 =23 y) 20) 
413 ~27 36 9 Oo] 
276 -19 7 ~12 08 
666 =37 35 - 2 “025 
TABLE 6. 


The largest values of diverrence in Table 6 are of the same order of 
marnitude as the findings of Baum [2] in his investigation of the relative 
sizes of the terms in the Baur and Phillip form of the tendency equation. 

The quantitative investigation of the present series of maps shows that 


the upver air reports are neither numorovs enoush nor accurate enouch for an 


(28 ) 
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evaluation of the terms t-VV or VV-t, on a quantitative basis. An error 
of orly one meter rer second in Vy can give an error of approximately 10 
millibars per 3 hours in the contribution of the term &-VYV . 
4e Synoptic Evidences of Convergence or Divergonoes 

Cenerally, it may be said that divergence is associatod with descend- 
ing air and resultant clear skies, while convergence implies ascending air 
and condensation formse Cyclonically curved and anticyclonically curved 
isobars aloft are frequently used to explain cloud forms and clear skies 
respectivelv. However, the mere fact of cyclonic or anticyclonic isobars 
does not imply either converrence or divergence above the friction layer, 
as can be seon by Figures l, 2, 5 and 6~¢ An absenco of all cloud forns 
implies either very dry air or widespread diversencese Similarly, wide- 
spread cloudiness implies widespread convergence, The types that the 
cloud forms assume depend on the convective stability of the layers con-« 
cerned in the lifting precess. Thus stable air gives stratiform clouds 
while unstable air gives cuniliform clouds on being lifted. Many pilots, 
the author being one of them, have often wondered why the clouds are fre- 
quently spaced in the atmosphere in layers or sheets. ‘tihen different levels 
of the atmosphere show different magnitudes of diverrtence or convergence, it 


follows that cloudiness will also exist in layers or sheetse 


(29) 
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